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Abstract 

Multi-wavelength studies of energetic solar flares with seismic emissions have 
I — , revealed interesting common features between them. We studied the first GOES 

X-class flare of the 24th solar cycle, as detected by the Solar Dynamics Observa- 
tory (SDO). For context, seismic activity from this flare (SOL2011-02-15T0f :55- 
X2.2, in NOAA AR 11158) has been reported in the literature (Kosovichev, 



C/2 



> 



m 



2011; Zharkov et ai, 2011[ ). Based on Dopplergram data from the Helioseismic 



. ' and Magnetic Imager (HMI), we applied standard methods of local helioseismol- 

"^ , ogy in order to identify the seismic sources in this event. RHESSI hard X-ray 

data are used to check the correlation between the location of the seismic sources 

O , and the particle precipitation sites in during the flare. Using HMI magnetogram 

■ ' data, the temporal profile of fluctuations in the photospheric line-of-sight mag- 



netic field is used to estimate the magnetic field change in the region where the 
seismic signal was observed. This leads to an estimate of the work done by the 
Lorentz-force transient on the photosphere of the source region. In this instance 
this is found to be a significant fraction of the acoustic energy in the attendant 
[~-^ . seismic emission, suggesting that Lorentz forces can contribute significantly to 

^^ ' the generation of sunquakes. However, there are regions in which the signature 

CJN , of the Lorentz-force is much stronger, but from which no significant acoustic 

emission emanates. 
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1. Introduction 

Since the discovery of a fiare-induced sunquake by IKosovichev and Zharkoval 
(|1998p many other works had been written about the detection of flare-generated 



seismic waves in the Sun (Zharkova and Zharkov, 2007 Donea et al., 2006 



Beshu-Ionescu et al, 2006 Moradi et al, 2007 Martinez- Ohveros et al, 2007 



Martinez-Ohveros, Moradi, and Donea, 2008). Sunquakes are defined as the 



observed expanding ripples on the solar photosphere induced by perturbations, 
driven by solar flares. Sunquakes carry energy and momentum into the solar 
photosphere and interior from the photospheric footpoints of a flare. They rep- 
resent acoustic waves that travel deep into the Sun and eventually refract back 
to the surface after a certain time. This is observed as a close-to-circular pattern 
of ripples in helioseismic Dopplergrams, expanding outward from the footpoints. 
The parts of the transient that travels along shallower ray paths come back to 
the surface nearer to the source, and sooner after the excitation, than that which 
penetrates deeper beneath the photosphere. The overall picture of this process 
is that of an outwardly expanding patterm of surface ripples seen up to about 
an hour after the impulsive phase of the flare. 

The importance of non-thermal particles (electrons and/or protons) interact- 
ing with the solar atmosphere has been already shown (Zharkova and Zharkov, 



2007; [Donea et al, 2006[ [Moradi et al, 2007[ |Martinez-01iveros et al, 2007D . 



The study of their acceleration, propagation and precipitation have proven to 
be a key area in flare seismicity research. A solar model divides the flare process 
into several stages: acceleration, injection, propagation, trapping, precipitation 
and energy losses. Each of these stages has a distinctive type of electromagnetic 
emission and timing associated with the physical processes occurring at that 
stage. Therefore, studying the flare emission at different wavelengths gives im- 
portant information about the physical conditions under which a sunquake can 
be generated. 

Recent studies of the acoustic emission of solar flares have opened a wide 
range of diagnostic possibilities and applications in helioseismology and in our 
understanding of these events ([Kosovichev and Zharkova, 1998 Donea, Braun, 



and Lindsey, 1999; Donea and Lindsey, 2005[ ). How sunquakes are produced is 



still poorly understood, but, several mechanisms explaining their generation have 
been suggested. Some of these hypotheses rely in the energy and/or momentum 
that electrons may transfer to the solar interior, directly (collisions) or indirectly 
(back- warming radiation). These ideas are based on the close relationship evident 
bet ween hard X-ray sources and the acoustic signatures ( Donea and Lindsey, 
2005; [Zharkova, 2008). Another hypothesis was given by Hudson, Fisher, and 



Welsch (2008), proposing that the magnetic field is itself the direct driver of 
acoustic activity during flares. In this case step-like changes of the magnetic field, 
as those observed by [Sudol and Harvey] (|2005p , do work on the photosphere via 



the Lorentz force. [Hudson, Fisher, and Welsch[(|2008p estimated the magnitude of 
this work and found it to be comparable to the values reported in the literature. 
In this work we perform a magneto-acoustic study of the first GOES X-class 
flare of the 24th solar cycle trying to assess if the magnetic field change during 
the flaring process could be an effective mechanism for the sunquake generation 
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or if is necessary to consider more than only one mechanism acting during the 
flare in order to produce its sunquake phenomena. 



2. Observations and Analysis 

On 15 February 2011 the active region designated NOAA AR11158 hosted a X2.2 
solar flare (SOL2011-02-15T01:55 X2.2) which began at 01:43:44 UT, reaching 
its maximum at 01:55:30 UT and ended at 02:27:32 UT, the flrst GOES X- 
class flare of the 24th solar cycle. To analyze this event we used Solar Dynamics 
Observatory/Helioscismic and Magnetic Imager (SDO/HMI) velocity and mag- 
netic data which consist of full-disk observations in the photospheric line Fe I 
6I73A, obtained with a cadence of 45 seconds ( jSchou et ai, 20f0| ). Aditionally, 
reconstructed images from the Reuven Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI), observed in the 12 — 40 and 40 — 100 keV energy bands, with 
a spatial resolution of ^ 1 arcsec, were used for the hard X-ray emission of the 

flare fsee llTurford et a/. I ([2002]) ). 

Seismic activity from this X2.2 flare was reported bv IKosovichevI (j2011|) . In 
this work, it is shown that the helioscismic waves were driven by a perturbation 
in the early impulsive phase, observed prior to the hard X-ray (50 — 100 keV) 
impulse, and probably associated with atmospheric heating by relatively low- 
energy electrons (~ 6 — 50 keV) and heat flux transport. However, low-energy 
electrons peak at the same time as high-energy electrons (this is the "soft-hard- 
soft" phenomenon), so there is no way to distinguish them. Heat flux transport 
is a strange concept in this situation in general, since the impulsive phase of 
a flare is so nonthermal and explosive. Aditionally, the seismic source is asso- 
ciated with one of the perturbations located along the flare ribbons and the 
anisotrop of the wavefront is explained in terms of the movement of this source 
(see IKosovichevI I201ip . Recently IZharkov et al.\ (|201ip analyzed the same event 
using acoustic h olography ([Donea, Braun, and Lindsey, 1999[ Donea, Lindsey, 
and Braun, 2000) showing the existence of a secondary weak seismic source. Also 
finding a spatial correlation between the seismic sources with the endpoints of 
a sigmoid, indicating the presence of a fluxropc in this two-ribbon flare. Based 
on this, IZharkov et al.\ (|201ip proposed a phcnomcnological scenario where the 
seismic sources are located at the footpoints of an erupting fluxropc. 

2.1. Helioscismic analysis 

Like IZharkov et al.\ (|201ip , our helioscismic analysis approach for this event is 
based on acoustic holography. Basically this technique corresponds to a phase- 
coherent reconstruction of acoustic waves observed at the solar surface into the 
solar interior to render "images" of sub-surface sources that have given rise to the 
surface disturbance. The main computations in holography are the "ingression" 
and "egression" power maps. These two quantities are estimates of the wave-field 
in the solar interior; the ingression is an assessment of the observed wave-field 
converging upon the focal point while the egression is an assessment of waves 
diverging from that point. 
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Using this computational technique, we reconstruct egression power maps of 
the flaring region over different ranges of the spectrum. Figure 1 shows an spec- 
trum of seismic emission histogram of the acoustic energy distribution derived 
from egression-power snapshots in the impulsive phase of the flare, because in 
this phase is when the sunquake mechanism is assumed to act. For the stronger 
source (source 1; see figure [D for source locations) most of the power is radiated 
in the 2.5-4.5-mHz spectrum; adding the 3- and 4-mHz energies gives us an 
acoustic energy of ^ (1.8 ± 0.15) x 10^^ erg. The fainter source (source 2) was 
visible in the 6 mHz band with an estimated acoustic energy of ^ 3 x 10^^ erg. 



Spcc;rum of Seismic Emission 



-From Source 1 into Pupi 
-From Quiet Sun, Total 
-From Quiet Sun into Pupil 



From Source 2 into Pupi 



Frequency(mHz) 

Figure 1. Energy spectrum of seismic emission from the flare of 15 February 2011 in 
NOAA AR11158. The stronger source (source 1) was visible in multiple frequencies, while 
the fainter source (source 2) was only clearly visible in the 6 mHz band. See figure[2]for source 
identification. 

The green curve, then, represents the egression power for the limited pupil. 
The blue curve is the total acoustic power. We integrate the blue and green 
curves and take the ratio as an estimate of the fraction of energy emitted by 
acoustic sources that returns to the surface in the limited pupil of the egression 
computation. These ought to be the same if all of the acoustic power emitted by 
the source went into the pupil and no other noise contributed to the egression- 
power signature. The reason the blue and green curves are as similar as can 
be seen in the Figure [T] is that most of the power from a source of the size of 
that represented by the egression-power signature in Figure [5] is captured by the 
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pupil — in spite of the fact that the rest of the Sun's surface offers a much greater 
arecjll- 




■-1.0 



50 Mm 



Figure 2. Acoustic source distribution for the flare of 15 February 2011 in NOAA AR11158. 
Upper panels show the 2.5-4.5-niHz and 5.5-6.5-niHz egression power maps normalized to 
unity in the mean quiet Sun. Center left frame shows the continuum image of the flaring 
region while the middle right frame shows a threshold mask whose value is unity (dark) in the 
source region and null (green) elsewhere. Lower left shows the pre flare line-of-sight magnetic 
field in kG. Lower right shows the base-ten logarithm of the mean square line-of-sight magnetic 
variation in the 2.5-4.5 mHz spectrum during the impulsive phase of the flare. Arrows labeled 
"SI" and "S2" show the locations of compact acoustic sources whose energy spectra are labeled 
"Source 1" and "Source 2," respectively, in Figure [T] The lack of a significant signature at the 
location labeled "S2" in the upper left frame signifies that excess seismic omission from Source 
2, while significant in the 5.5-6.5-mHz band, is insignificant in the 2.5— 4.5-mHz band. 



TThis is a result of strong refraction beneath the solar surface due to a sound speed that 
increases rapidly with depth (see [Lindsey and Braun| II2000II ) . If the sound speed were uniform 
in the solar interior, nearly all of an acoustic power emitted by the source would miss a small 
pupil such as that applied in the egression-power computations. 
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Figure 2 shows a snapshot of the 2.5-4.5 mHz egression power (upper left) during 
the impulsive phase of the flare cospatial with the line-of-sight magnetic field 
(lower left) and the base-ten logarithm of the mean square magnetic variation in 
the 2.5-4.5 mHz spectrum (lower right). The middle right frame presents a mask 
that identifies what we regard to be the main source region, obtained simply by 
applying a threshold to the egression-power signature of the acoustic emission 
evident in the upper left frame. 

2.1.1. Magnetic field variation 

According to IFisher et al.\ (I2012[) , given a change in the photospheric magnetic 
field, we can estimate the components of the Lorentz force resulting from mag- 
netic restructuring. The temporal change of the force per unit area is given 

by 

^^^ ' ^ dAl- (Bl Bl) , (1) 



dt Stt J a dt 

where B±^ and B, are the two components of the magnetic field, and the integra- 
tion is made over the area A where the flare-driven fleld changes occur, giving 
the total Lorentz force on the uppermost layer of the solar envelope. We can 
rearrange the last equation in the form 

Following the argument given by Fisher et. al (2011), we can assume that the 
flclds are observed to change over a time interval Ai, and perform a temporal 
integration in the last expression in order to find the changes in the component 
of the Lorentz force, AF^: 

AF, = -^ / dA (Bi • ABi - B.AB,) . (3) 

As we are interested in the magnetic perturbation over the photospheric surface, 
we consider the inward force resulting from the magnetic restructuring during 
the flaring process: 

^F^.in^ard = -AF, = ^ f dA {B,AB, - Bi • ABx) . (4) 

In order to quantify the Lorentz-force contribution to the energetics of the 
sunquake, some theoretical asumptions must be made taking into account the 
available data for the estimation. First of all we will assume that ABz is null, 
which it would have to be in the case of magnetic flux frozen into a photospheric 
medium that is rigid to the extent that we cannot see a change in the geometry 
of the magnetic region over the acoustic relaxation time, Tac ^ c/2H '^ 40 s. 
Then B^AB^ — B^ ■ ABj^ is simply — B^ • ABj^, which we will suppose is 
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of the same order of magnitude as SiosABios- This can possibly be mistaken. 
For example, AB could have a magnitude of five times that of AB\os if AB^ 
were within about a fifth of a radian of being perpendicular to the line of sight. 
The probability of this would appear to be inversely proportional to the angle 
to within which AB would have to be perpendicular to the line of sight. 

What follows, then, is an estimate of B\osAB\os, but the designation "(los)" 
is left off from this point on. We compute the mean value of the line-of-sight 
magnetic field over the mask shown in the Figure [21 to find that 



{B) 



500 G, 



(5) 



where the angular brackets indicate the spatial average of their contents over 
the region identified by the source mask (sec Figure [21 top right panel) . 
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Figure 3. Dopplcr (upper panel) and line-of-sight magnetic field (lower) integrated over the 
source region defined by the threshold mask represented in the upper left panel of Figure[2]are 
plotted as a function of time, t. Vertical green linos at time zero identify the impulsive phase 
of the fiare. Red fang-like profile in the upper frame represents the Doppler profile of the fiare 
uncontaminated by the 5- minute oscillations. Dashed red lines in the lower frame identify what 
we propose to recognize as a step-function transient component of the line-of-sight magnetic 
profile. 



SOLA: SP_Alvarado_et_al_2011.tex; 27 February 2013; 17:28; p. 7 



Alvarado-Gomez et al. 



Variations in the Doppler and line-of-siglit-magnetic signatures spatially inte- 
grated over the source region are plotted in Figure [3] The top frame shows 
the Doppler signature (red shift) integrated over the source region. This shows 
the regular 5-minute oscillations and a significant transient behavior during the 
impulsive phase of the flare, which occurs at time zero. The bottom frame shows 
the line-of-sight magnetic field strength, B, integrated over the same region 
minus the 500 G mean. It also shows the transient behavior of the field. We 
suggest that this profile, in the time frame of the impulsive phase, has an initial 
excursion of some 37 G that is mostly ephemeral; the excursion returns to a 
value that is ~12 G above from the pre-fiare trend, which was increasing at a 
rate of 0.4 G/min, indicated by the slope of the dashed red lines in the bottom 
frame of Figure [31 

Now, we take the significant magnetic transient to be the stcp-function-like 
component shown by the displacement between the dashed red lines in the lower 
plot, the displacement between these being 12 G: 



A(B) = 12 G, (6) 

which is a small change compared to previous magnetic field restructuring stud- 
ies during energetic flares; for instance, [Sudol and Harvey] (|2005p report, for a 
sample of 15 X-class solar flares, a mean value of 90 G for the magnetic field 
change, using Global Oscillation Network Group (GONG) magnctograms, while 
IZharkova et a/.l ( |2005p . using Michelson Doppler Imager (MDI) magnctograms, 
report an absolute change of 120 G for a GOES X4.8 class fiare. 

In principle, these changes are in line with a model of magnetic-energy release 
that requires a "permanent change" in the magnetic configuration to drive a 
sunquake. Physically speaking, the energy released by a coronal magnetic field 
leaves the coronal field permanently changed in some respect, so that B'^/Sir in- 
tegrated over the corona reflects the reduced magnetic energy following the fiare 



( Hudson, 2000 1 . The magnetic virial theorem requires a significant manifestation 
of this changed configuration in any surface, S, that effectively "surrounds" the 
region of reduced magnetic energy, so long as the field is force free on S and 
throughout its interior. This strongly suggests that such a change will have a 
manifestation in the photospheric magnetic field, through which the magnetic 
flux passes. Because of this, a step-function component in the photospheric 
magnetic signature marking the impulsive phase of the fiare is of particular 
interest as a prospective source of work done on the solar interior acoustic field. 
An ephemeral change that is reversed once the impulsive phase of the flare has 
run its course is suspect in that it might be no more than the signature of a 
perturbed radiative environment during the impulsive phase. Of course, it is 
possible that the radiative environment itself is permanently changed after the 
impulsive phase of the flare. Moreover, the required change need not occur in 
the neighborhood of the acoustic source to satisfy the magnetic virial theorem. 
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In accordance with the foregoing considerations, we take the transient mag- 
netic force density (per unit area, hence a pressure) to bql| 

Apm = ^{BAB) ~ (B)A(B) = 480 ^. (7) 

47r cm'' 

Counting the number of unity-value pixels in the mask, we find an area of 

A = (105 pix) dx dy, (8) 

with pixel dimensions 

dx ^ dy = 5x IQ-^Rq = 0.35 Mm, (9) 

hence, 

A = 12.8 Mm^ = 1.3 x 10^^ cm^ (10) 

This leads to a total force transient integrated over the source region of 

AF^ = Ap^A = 3.12 X 10^^ dyne. (11) 

Now, to estimate the work done on the acoustic field by this force, we need 
to determine how far, Az, it depresses the photosphere in the time frame of 
the acoustic relaxation time, roughly 40 s in the quiet photosphere. One way 
to do it is to integrate over the fang-like profile plotted in red in the top panel 
of Figure [Sj supposing that, and admitting significant errors, the displacement 
uncontaminated by 5-minute oscillations would look something like this. We 
take the area thus subtended to be 



1 
Az = -basexhcight = 0.5x225 sx320 m/s = 36 km = 3.6x10^ cm, (12) 

where the base and height of the profile represent the duration (225 s) and 
velocity (320 m/s), respectively, of the photospheric depression. We now take 
the work done by the force over the displacement supposing that the mean 
imbalance of the force as the displacement progresses is half the initial transient, 
AF, above; then, 

AW = - AF^ Az = 1.12 X lO^Srg. (13) 

The energy estimate expressed in equation ((13)) is less than the ~ (1.8 ± 0.15) x 
10^^ erg based on the helioseismic signature of the flare (source 1). Given the 



t Because -Bios varies considerably across the source region, the mean of BAB can be consider- 
ably different from the product of the means used as a proxy in equation lO . The error could 
be of order 100%. However, we think it is unlikely to be anything like an order of magnitude. 
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other relatively rapid changes in the line-of-sight magnetic signature evident in 
the Figure [31 the possibility ought to be considered that the force transient is 
greater than the step-function component. If we suppose that it is something like 
the full 37 G range of the excursion seen between in the time interval —1.0 to 
+2.5 min, for example, this would increase the estimate of AW proportionately, 
to 3.36 X 10^^ erg. This value reaches ^ 18% of the egression-power estimate for 
the analyzed source, representing too considerable a fraction of it to be regarded 
as negligible; it must be recalled at this point that this estimation is considering 
only the line-of-sight component of the magnetic field, and the total energetic 
contribution from the magnetic field restructuring can be even higher than the 
value presented here. 

2.1.2. Hard X-ray emission 

Using RHESSI hard X-ray data we determined the spatial location of the hard X- 
ray footpoints in four different time intervals, as is shown in the Figure 21 during 
the impulsive phase where the transients in both, the magnetic and Doppler, were 
observed. These observations revealed a complex structure in the distribution 
of footpoints during this phase, suggesting multiple magnetic loop formation 
and movement along the flare ribbons. This process was likely to end with the 
relaxation of the magnetic field lines and the precipitation of particles in the 
endpoints of one the flare ribbons. The left white arrow indicates the center of 
the analyzed region (mask in the Figure [2]). 



Doppler Difference 01:49:57 



Running Difference 02:08:42 





180 200 220 240 

X (orcsecs) 



Figure 4. HMI dopplcrgram differences at the beginning of tlie impulsive pliasc (left) and at 
the sunquake wavefront detection time (right), combined with RHESSI CLEAN contour plots, 
with levels 50%, 60%, 70%, 80% and 90% of the maximum detected intensity, in the 40 - 100 
and 12 — 40 keV bands, integrated during the intervals shown. The white arrows indicate the 
location of the endpoints of one of the flare ribbons. The left arrow points to the center of the 
analyzed region. 



SOLA: SP_Alvarado_et_al_2011.tex; 27 February 2013; 17:28; p. 10 



Magneto-Acoustic Energetics Study 



3. Discussion 

We have studied the first GOES X-class flare of tlie 24tli solar cycle detected by 
HMI /SDO, the 15 February 2011 X2.2 flare posovichev, 2011[ Zharkov et ai, 
2011), confirming the two induced sunquakes by this event. Hudson, Fisher, and 
Welsch ^2008p and lFisher et a/.l(|2012p prescribe a means to estimate the contri- 
bution of the Lorentz-force transient to such a sunquake based upon transient 
vector magnetic signatures across the impulsive phase of the flare. What we have 
undertaken here is a preliminary, order-of-magnitudc estimate of the Lorentz- 
force contribution based on just the line-of-sight magnetic signature, pending 
the availability of reliable vector magnetograms from HMI. The full contribution 
based on the vector-magnetic measurements could be greater or less than this 
estimate, but basic statistics favor an expectation that is at least somewhat 
greater. Hence, if B and AB could be regarded as isotropically distributed, then 
the rough estimate we derive here is most likely somewhat of an underestimate. 
Our estimate of the work done by the Lorcntz force contribution based on 
just the line-of-sight transient magnetic signature is (1.12 ± 0.1) x 10^^ erg, 
approximately 6% of the 1.8 x 10^^ erg we estimate for the total energy in the 
acoustic transient. This alone is a significant fraction, and if we knew it to be 
no more than this, it would be dangerous to regard it as negligible for general 
purposes. In fact, if we revive the proposition that the vertical component of 
AB_L, is null, then it should follow that 

ABios = pAB • 6, (14) 

where AB is the componet of B^ = B that lies in the local vertical plane 
containing the line-of-sight direction, 6, and p is the sine of the angle between 
6 and the vertical direction, z. For AR11158 at the time of the flare, p was 
0.267. Hence, under the assumption that AB is horizontal, it follows that the 
real value of the component of AB to which the line of sight is sensitive is 
AB/p = (12/0.267) G = 45 G. Replacing the 12 gauss estimate of equation (|13p 
by 45 gauss increases the 10^^ erg estimate to 3.75 x 10^^ erg. This is about 23% 
of the 1.8 X 10^^ erg acoustic energy we estimated. This remains more likely an 
underestimate than an overestimate, and exceeds the acoustic energy in the (5-7 
mHz) where the mapping of the egression-power signature is most conspicuous 
(as it has been in previously detected sunquakes, see |Donea and Lindsey([^005p . 
At the same time, maps of the mean square line-of-sight magnetic transient 
(see lower right panel in Figure[l]for this in the 2.5-4.5 mHz spectrum) show rms 
variation in AiJios up to nearly 100 G with no correction for p in regions far from 
Sources 1 and 2 and from which no significant acoustic emission emanates. More- 
over, RHESSI, the hard X-ray emission of the event show a complex structure of 
footpoints with multiple locations of particle precipitation during the impulsive 
phase, without a clear spatial correlation between these precipitation sites and 
the location of the acoustic sources. Whether or not Lorentz-force transients 
contribute significantly to flare acoustic emission, then, appears to remain an 
open question, and this seems to involve some very interesting physics that has 
yet to be identified. 
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These results also show the need for a more detailed study that utilizes high- 
resolution vector magnetic observations and includes the acoustics of heating 
of the outer atmosphere by accelerated particles. Crucial to such a study are 
realistic numerical simulations of the physics of prospective contributors to flare 
acoustic emission. 
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